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Abstract 

We consider the use of polarization properties as a means to discriminate between 
Synchrotron Self-Absorption (SSA) and Free-Free Absorption (FFA) in GHz-Peaked 
Spectrum (GPS) sources. The polarization position angle (PA) of synchrotron radia- 
tion at high frequencies for the optically thin regime is perpendicular to the magnetic 
field, whereas it is parallel to the magnetic field at low frequencies for the optically 
thick regime. Therefore, SSA produces a change in PA of 90° across the spectral 
peak, while FFA does not result in such a change. We analyzed polarization data 
from VLA observations for six GPS sources to see if such a change in PA was present. 
Our results indicate that there is no significant evidence for 90° change in PA across 
the spectral peak, suggesting that FFA is more likely than SSA for low-frequency 
cutoffs in these sources. 
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Introduction 

GHz-Peaked Spectrum (GPS) sources are characterized by a simple convex radio spec- 
trum with a peak at GHz frequencies. The common properties of bright samples of GPS 
sources are: small size (< 1 kpc), powerful radio luminosity (L > 10 45 erg s _1 ), and ap- 
parently low variability (e.g., 000 [cite]cite.Stanghellinil998Stanghellini et al. (1998); 000 
[cite]cite.O'Deal9980'Dea (1998)). GPS sources associated with quasars tend to show a 
complex morphology in VLBI images, while those associated with galaxies are likely to be 
double or triple radio sources (000 [cite]cite.Phillip-Mutell980Phillips, Mutel, (1980); 000 
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[cite]cite.O'Deal9910'Dea et al. (1991)). 

The spectral shape of GPS sources is distinct from that of extended radio sources. 
Generally, extended radio galaxies have steep spectral profiles, in contrast to the flat spectral 
profiles of quasars. GPS sources have steep spectra at high frequencies, as do radio galaxies, 
which indicates optically thin synchrotron radiation. However, the low-frequency cutoff in the 
spectral profile is peculiar to GPS sources. There has been a considerable debate regarding the 
interpretation of this low-frequency cutoff. 

The low-frequency cutoff is ascribed by some to Synchrotron Self- Absorption (SSA). 
SSA is usually discussed in terms of the equipartition brightness temperature, T B (assuming 
equipartition of energy between the radiating particles and the magnetic field), which is derived 
from the flux density and frequency of the spectral peak, and the optically thin spectral index. 
|e]cite.Readheadl996Readhead et al. (1996) ([cite]cite.Readheadl996Readhead et al. (1996)) 
showed that the observed peak brightness temperature is consistent with the equipartition 
brightness temperature for the lobes of 2352+495. 000 [cite]cite.Snellen2000Snellen et al. (2000) 
suggested SSA based on correlations between the peak frequency and the angular size, and 
between the peak flux density and the angular size. 

On the other hand, ||e]cite.Bicknelll997Bicknell et al. 

(1997) ([cite]cite.Bicknelll997Bicknell et al. (1997)) suggested that the observed anti- 
correlation between z/ m and size can be explained by Free-Free Absorption (FFA). VSOP and 
multi-frequency ground VLBI observations of OQ 208 have been used to infer that FFA is the 
more likely mechanism, with the two radio lobes surrounded by a cold dense plasma which 
causes thermal FFA (000 [cite] cite. Kameno2000Kameno et al. (2000)). Assuming intrinsically 
symmetric double lobes, the apparent asymmetry of opacities towards lobes was explained as 
being caused by the differences in path lengths through the plasma, which has an estimated 
electron density and electron temperature of 600 < n c < 7 x 10 5 and 10 4 < T c < 6 x 10 7 , 
respectively. In addition, FFA is favored for other GPS sources, such as 0108+388 (000 
[cite]cite.Marr2001Marr et al. (2001)) and NGC 1052 (000 [cite]cite.Kameno2001Kameno et 
al. (2001)). 

So far, the origin of the low-frequency cutoff has been discussed in terms of the spectral 
properties and brightness temperatures. However, since both FFA and SSA can produce similar 
spectral shapes, it is difficult to differentiate between the two mechanisms in this way. Here, we 
examine the possibility of discriminating FFA from SSA by using the polarization properties of 
sources. In section 2 we review the relevant polarization properties, and we describe the data 
analysis in section 3, followed by section 4 where the results for the sources studied to date are 
given. 

2. A New Test of Absorption Mechanisms Using Polarization Properties 



te]cite.Allerl970Aller (1970) ([cite] cite. Allerl970Aller (1970)) pointed out that the po- 
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larization position angle (PA) is perpendicular to the magnetic field when synchrotron emission 
is optically thin in terms of SSA, while PA is parallel to the magnetic field in optically thick 
regimes. This is explained as follows. The ratio of the specific intensities, Ib/Ia, determines the 
PA, where I\ and I& are specific intensities parallel and perpendicular to the magnetic field, 
respectively. The specific intensity is approximately proportional to the emissivity, e, when the 
emission is optically thin, and can be described as 
J b^b_ 3 7 + 5 

I a e A 2 K J 

Here, we assume that the number density of relativistic electrons follows a power law, N(E)dE = 
KE~^ dE, as a function of energy E where K is a constant. When the emission is optically 
thick, the specific intensity is also inversely proportional to the absorption coefficient, k, so that 
h e B /K B 37 + 5 



< 1- (2) 
I a e K /n A 3 7 + 8 V ' 

Comparing equations (1) and (2), we can see that PA changes by 90° between optically thick and 

optically thin regimes. The optical depth of synchrotron radiation decreases with increasing 

frequency, so that the radiation will be optically thick at low frequencies and optically thin 

at high frequencies. The spectral peak occurs at the frequency where the optical depth is 

approximately unity. Therefore, the 90° jump in PA appears across u m . 

Such a 90° jump of PA is found in some BL Lac objects. Absorption in BL Lac ob- 
jects is thought to be due to SSA (Q)0 [cite] cite. Readheadl996Readhead et al. (1996); 000 
[cite] cite. Shenl999Shen et al. (1999)). Figure 1 shows (a) the total flux densities and (b) the 
PAs of the BL Lac object 0212+735 over a period of 4 yr from observations at the University of 
Michigan Radio Astronomy Observatory (UMRAO). The PA at 14.5 GHz jumped by 90° after 
MJD 40100, while at 4.8 GHz it remained almost constant. The total flux densities indicate 
that the emission was optically thick between 4.8 and 14.5 GHz before MJD 40100, but became 
optically thin at 14.5 GHz. This variability is thought to be due to the expansion of a jet com- 
ponent. This synchronous change of the opacity and PA suggests a causal relationship between 
them, implying SSA rather than FFA in BL Lac objects. Another example of a 90° jump is 
reported for the BL Lac object OJ 287 (000 [cite]cite.Gabuzda2001Gabuzda, G'omez (2001)). 
They detected a 90° jump in PA for the core between 5 GHz and 22 GHz, and interpreted it 
as the transition from the optically thick to optically thin domain. 

The above properties of PA can be used as a test of the absorption mechanism in GPS 
sources. If the low-frequency cutoff is caused by SSA, we expect to observe a 90° jump in 
PAs across the spectral peak frequency. Otherwise, in the case of FFA, no such 90° jump will 
be seen. (This, of course, assumes that the emission from optically thick regime dominates 
the polarized emission, an assumption which is easily tested by observations with sufficient 
angular resolution to resolve the contributions to the polarized flux from the various source 
components.) Based on polarization data just below z/ m , ||e]cite.Stanghellinil998Stanghellini et 
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al. (1998) ([cite] cite. Stanghellinil998Stanghellini et al. (1998)) mentioned the possibility that 
the polarized emission may not come from the region responsible for the optically thick emission, 
or that the turnover is not due to SSA. Because GPS sources are generally not time variable, 
we combined multi-frequency observations around v m from various epochs to investigate the 
existence of the 90° jump. 

3. Data Analysis 

3.1. The Samples from VLA Observations 

We analyzed the PAs of six GPS sources observed with the VLA (000 
[cite] cite. Saikial998Saikia et al. (1998); 000 [cite]cite.Stanghellinil998Stanghellini et al. 
(1998)). Our selection criteria were: 

(1) to have more than 2 polarization observations in both optically thick and thin fre- 
quency regimes, 

(2) to have a percentage polarization higher than 0.1% at 5 GHz. 

The sources 0248+430, 0711+356, 0738+313, 0743-006, 1143-245 and 2134+004 met these 
criteria. They have all been identified as quasars, and their properties are summarized in 
table 1. These sources have between 4 and 6 polarization observations at different frequencies. 
The PAs and percentage polarizations are given in table 2, and the PAs of the six sources as 
a function of squared wavelength (A 2 ) are plotted in figure 2. The vertical lines in these plots 
indicate u m . 

3.2. Statistical Study 

We analyzed the distribution of the PAs statistically. In the case of FFA, PA follows the 
Faraday Rotation formula across the spectral peak, which can be fitted by 



at both optically thick and thin frequencies. Here, 6^ is the PA observed at A^, and A m is 
wavelength corresponding to u m . The free parameters RM and 8q are the Faraday Rotation 
Measure and the intrinsic PA, respectively. In the case of SSA, PA is fitted as follows: 



The sign "±" indicates jump in PA between the optically thick and thin regimes. Figure 3 
shows an example of PAs for 1143—245, with and without ±90° offsets for the data points 
above the peak frequency. We fit both the FFA and SSA models to the data and evaluated the 
residuals for each source. Here, we introduce " dist(9) n as the distance of PA from the fit, 



< k > FFA = 9 + RMXl 



(3) 




(4) 




(5) 
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F(,ssa,^pa) = (9) 



= sin 2 0. (6) 

This function has a 180° period, and is thus free from the 180° ambiguity in PA. When 9 is 
sufficiently small, dist(9k— < 9k >) is approximately equal to (9k— < 9k >) 2 . Here, we define 
the residual, S, as 

s _ dist(9 k - < 9 k >) ^ 
k a k 
We use the ratio, F rat , defined as 

P - SssA 

>->FFA 

to determine which absorption mechanism is more likely. Unless there is any significant differ- 
ence of likelihood between the FFA and SSA models, as a first-order hypothesis, the ratio, F rat , 
would follow the F-distribution, 

Xssa/^ssa 
Xffa/ ^ffa 

where i^ssa an d ^ffa are the numbers of degrees of freedom corresponding to the \ 2 values, 
Xssa an d Xffa; respectively The ratio of variances, 

^ssa _ ^SSA 

xl FFA SWa 

is also distributed as F rat . Evaluating F rat , one can thus discriminate between two models. 
Simply, if F rat > 1, then FFA is more likely. 

4. Results 

In table 3 and figure 4, we show the results of statistical analysis. The value of F Tat 
is greater than 1 for all sources except 2134+004. This implies that FFA is more likely than 
SSA. For 0711+356, F rat = 8.586 indicates that the residual of the FFA model is significantly 
less than that of the SSA model. This means that a simple SSA model cannot be accepted as 
the cause of the spectral cutoff below 1.6 GHz in this source. In the cases of 0738+313 and 
0248+430, F rat = 2.431 and F rat = 2.067, respectively, again indicating that the FFA model 
is more likely. For 1143—245, _F rat = 1.304. This is because of the high RM in this source, 
which results in small values of S for both models (see figure 3). For 0743—006 and 2134+004, 
because we have only four data points, the degree of freedom is very small. Consequently, the 
values of S for both SSA and FFA are small, and F rat is close to 1, as can be seen in figure 4b. 

5. Discussion 

Although five of the six sources have _F rat > 1, the case for FFA is not compelling, in part 
because of lack of observed frequencies, particularly around z/ m . The SSA model predicts not 
only the 90° jump in PA, but also the change in the fractional polarization as a function of the 
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wavelength. In the optically thin regime, the degree of polarization decreases with decreasing 
frequency, and becomes null at v m (000 [cite] cite. Allerl970Aller (1970)). It increases again at 
the optically thick regime, but the polarization degree is ~ 1/10 of that at the optically thin 
regime. The percentage polarization is plotted as a function of A 2 for the six sources studied 
here in figure 5. The variation in the degree of polarization expected from the SSA model 
cannot be seen for any source except 0738+313. This also implies that SSA is unlikely to be 
the cause of the low-frequency cutoff. 

For 2134+004, F rat <1 results from the low spatial resolution against the complex struc- 
ture of the source components. In fact, 2134+004 has two components, which have different 
RMs (000 [cite]cite.Taylor2000Taylor (2000)). In such a case, our method should be applied to 
the individual components. The RM asymmetry of 2134+004 is easily explained by means of 
the FFA model. Suppose that two intrinsically identical lobes are surrounded by dense plasma 
to produce FFA. Because there is usually a difference in the path lengths to the lobes in the 
plasma, the RM differences reflect the differences in the path length. Such a model allows us 
to estimate the electron density and the magnetic field (000 [cite]cite.Kameno2000Kameno et 
al. (2000)). 

In conclusion, studies of the polarization properties of six published GPS source have 
yielded no convincing evidence for the 90° change in PA across the spectral peak expected from 
SSA models. This and the distribution of fractional polarization as a function of wavelength 
lead us to conclude that FFA may be the dominant cause for the low-frequency cutoffs in most 
of these sources. 



This research made use of data from the University of Michigan Radio Astronomy 
Observatory, which is supported by funds from the University of Michigan. We thank Philip 
G. Edwards for much advice. 
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Table 1. Sample of GPS sources. 



Source 


z 


m 


S^ghz (Jy) 


fmax (GHz) 


S'max (Jy) 


0248+430 


1.316 


15.5V 


1.24 


5.2 


1.27 


0711+356 


1.620 


17.0V 


0.861 


1.6 


1.42 


0738+313 


0.631 


16.1V 


3.66 


5.1 


3.82 


0743-006 


0.994 


17.5V 


2.05 


6.0 


2.12 


1143-245 


1.950 


18.5V 


1.40 


2.0 


1.69 


2134-004 


1.936 


16.8V 


8.59 


5.2 


8.59 



The columns are: source name, red shift, optical magnitude and filter/color, flux density at 5 GHz, observed 
peak frequency and peak flux density (Stanghcllini ct al. 1998). 
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Table 2. Polarization of GPS samples 



frequency 


1335 




1665 




4535 




4985 




8085 




8465 




source 


% 


PA 


% 


PA 


% 


PA 


% 


PA 


% 


PA 


% 


PA 


0248+430 


1.7 


+52 


1.5 


+23 


2.0 


-18 


2.0 


-27 


0.5 


-40 


0.4 


-44 


0743-006 


0.4 


-76 


0.5 


-7.0 


0.2 




0.2 




0.8 


+1.0 


1.0 


-18 


1143-245 


1.9 


+72 


2.1 


+80 


0.7 


-18 


0.7 


-27 


1.3 


-43 


1.4 


-46 


2134—004 


0.1 




0.1 




0.9 


+60 


0.9 


+43 


0.5 


-4.0 


0.5 


+0.0 


frequency 


1380 




1630 




4815 




4865 




8435 




8485 




source 


% 


PA 


% 


PA 


% 


PA 


% 


PA 


% 


PA 


% 


PA 


0711+356 


1.3 


+5.0 


1.0 


-33 


1.8 


+83 


1.8 


+80 


1.9 


+76 


2.0 


+79 


frequency 


1465 




1665 




4835 




14915 




14965 








source 


% 


PA 


% 


PA 


% 


PA 


% 


PA 


% 


PA 






0738+313 


0.4 


+30 


0.3 


+24 


6.3 


+45 


6.6 


+26 


6.7 


+25 







Note. Columus are source name, fractional polarization in %, PA (polarization position angle) in degree. 
The data were taken from Stanghcllini ct al. (1998), except for those of 0738+313 which were adapted from 
Saikia et al. 1998 



Table 3. F-test of the model fit 



Source 


name 


Degree of freedom 


F rat [=S(SSA)/S(FPA)] 


Significance level [%] 


0248H 


-430 


4 


2.067 (0.137/0.066) 


24.97 


071H 


-356 


4 


8.586 (0.103/0.012) 


3.04 


0738H 


-313 


3 


2.431 (0.073/0.030) 


24.24 


0743- 


-006 


2 


1.304 (0.043/0.035) 


45.31 


1143- 


-245 


4 


1.207 (0.038/0.029) 


40.15 


2134- 


-004 


2 


0.983 (0.046/0.047) 


50.43 



The columns are: source name, degree of freedom, i^ ra t , and significance level for F va t [% 
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Fig. 1. (a) Flux density as a function of MJD. The total flux density indicates that the emission had 
been optically thick from 4.8 to 14.5 GHz before MJD 40100, and then became optically thin at 14.5 GHz 
thereafter, (b) PA variation as a function of MJD. There is a 90° jump in PA at 14.5 GHz after MJD 40100. 
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Fig. 2. PA as a function of A 2 for the six GPS sources. There are four to six data points for each source. 
The vertical line shows the spectral peak wavelength. In most cases it seems that a straight line can be 
fitted throughout the peak wavelength, in support of the FFA model. 



1143-245 




Fig. 3. PA as a function of A 2 for 1143—245. Two data points at longer wavelengths arc also plotted with 
offsets of ±90°. Due to the high RM, these offset data could be fitted well by straight lines. 
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Fig. 4. (a) The vertical bar indicates F rat [=S(SSA)/S(FFA)]. These are all greater than 1, except 
2134+004. (b) Plot of S(SSA) against S(FFA). The dotted line indicates S(FFA)=S(SSA). For all sources, 
except 2134+004, the data points are above the dotted line. 
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Fig. 5. Fractional polarization as a function of A 2 . The vertical line shows the spectral peak wavelength. 
All except 0738+313 show essentially the same fractional polarization in both optically thin and thick 
frequencies, again in favor of FFA. 



